Localization of Dirac electrons by Moire patterns in graphene bilayers 
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We study the electronic structure of two Dirac electron gazes coupled by a periodic Hamiltonian 
such as it appears in rotated graphene bilayers. Ab initio and tight- binding approaches are combined 
and show that the spatially periodic coupling between the two Dirac electron gazes can renormalize 
strongly their velocity. We investigate in particular small angles of rotation and show that the 
velocity tends to zero in this limit. The localization is confirmed by an analysis of the eigenstates 
which are localized essentially in the AA zones of the Moire patterns. 
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Introduction - Graphene is a two dimensional carbon 
material which takes the form of a planar lattice of sp 2 
bonded atoms. Its honeycomb lattice consists in two sub- 
lattices which gives a specific property to the wave func- 
tion, the so-called chirality. The linear dispersion rela- 
tion close to the charge neutrality point implies that, in 
this energy range, the electrons obey an effective mass- 
less Dirac equation [TJ [5] [3]. The properties of elec- 
trons in graphene, deriving from the Dirac equation, 
are fundamentally different from those deriving from the 
Schrodingcr equation. In particular the quantum Hall 
effect is quantized with integer plus half values [U [5] and 
can even be observed at room temperature [5] . Another 
example of the unique behavior of Dirac electrons is the 
so-called Klein paradox which is intimately related to the 
chirality of their wavefunction. The Klein paradox is the 
fact that, in a one-dimensional configuration, a poten- 
tial barrier is perfectly transparent for electrons. As a 
consequence it is difficult to localize Dirac electrons with 
an electrostatic potential, although it can be of great in- 
terest to realize this confinement, in particular for the 
production of elementary devices [TJ, [5] • 

Therefore it is necessary to improve our knowledge of 
the behavior of Dirac electrons in various types of po- 
tentials. For disordered potentials specific behavior re- 
lated to the chirality of Dirac electrons, such as weak 
anti-localization effects, have been predicted theoreti- 
cally [9] and observed in epitaxial graphene [10), [TTJ [12] 
H5] . Recent STM measurements, performed on epitaxial 
graphene, have also confirmed the absence of backscat- 
tering related to the Klein paradox [T3] . For periodic 
systems several studies have been published concerning 
either a simple graphene sheet in an applied periodic po- 
tential or periodic bilayers [TSj |TB]. The bilayer system 
[Tfl IT51 Q1|J [50] has shown an unexpected and rich be- 
haviour. The well-known Bernal AB stacking leads to 
massive quasiparticules with quadratic dispersion close 
to the Dirac point. Yet it has been recognized recently in 
the scientific community that another important case is 
that of rotated graphene bilayers, which are observed for 



example in epitaxial graphene on the C-terminated face 
of SiC [17 . In a rotated bilayer the superposition of the 
two rotated honeycomb lattices generates a Moire pat- 
tern with a longer period. The two Dirac electron gases 
are then coupled by a periodic interaction, with a large 
supercell, which can restore a Dirac-like linear disper- 
sion. However fundamental issues are controversial. In 
particular Lopez dos Santos et al. |18j have shown, with 
a perturbative treatment of the interplane coupling, that 
the velocity can be renormalized and tend to decrease in a 
rotated bilayer with respect to its value in pure graphene. 
In contrary Shallcross et al. 19J conclude from ab initio 
calculations and from some general arguments that the 
velocity is unchanged in rotated graphene bilayers. 

In this letter we combine ab initio and tight-binding 
(TB) approaches and conclude indeed that there is a 
renormalization of the velocity for bilayers as compared 
to pure graphene. Most importantly the velocity tends to 
zero in the limit of small twist angles which means that 
this type of coupling is able to confine the two Dirac elec- 
tron gazes. We also show that the electronic wavefunc- 
tion tends to localize in regions of the bilayer which are 
locally similar to the AA arrangement. This new local- 
ization regime should be observable since angles as small 
as a fraction of a degree, for which we predict strong lo- 
calization effects, are found in some rotated multilayers. 
In particular large Moire patterns have been observed on 
STM image of graphene multilayers on top of SiC C-face 
[51] and also on graphite [25]. Our work demonstrates 
also that the perturbative theory of Lopez dos Santos et 
al. [18] is correct for angles greater than typically 3 de- 
grees but cannot deal with the remarkable localization 
phenomenon that occurs in the limit of small twist an- 
gles. 

Moire patterns and geometry of rotated bilayer.- Moire 
Patterns can be obtained in two cases: when two lattices 
with slightly different parameters are superimposed or 
when two identical lattices are rotated by an angle 9 [23] . 
The present situation corresponds to the later case. A 
commensurate structure can be defined if the rotation 
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FIG. 1: (color on line) Commensurate bilayer cell (n, m) = 
(6, 7) for a rotation of 6 = 5.08°. Full (dashed) line circle AB 
(AA) region. 

changes a lattice vector V (to, n) to V (n, m) with n, 
m the coordinates with respect to the basis vectors a\ 
(73/2,-1/2) and a 2 (y/3/2,1/2). The rotation angle is 
then defined as follows: 

n 2 + Anm + m 2 

cos9=—- — ■ (1) 

2{n z + lira + mr) 

and the commensurate cell vectors correspond to : 

t = V = ric?! + ma 2 ; t' = —ma\ + (n + m)a 2 . (2) 

The commensurate unit cell contains N = 4(n 2 + nm + 
to 2 ) atoms. It is important to notice that 9^0 results 
in a large cell -large n and to- and small \m — n\. Large 
commensurate cells can be also obtained for large angles 
9 ~ 30° then \m-n\ is large [23]. 9 = 60° is the perfect 
AB stacking and 9 close to 60° is obtained for n = 1 
(to = 1) and large to (n). 

For small angle, if the rotation axis, perpendicular to 
the planes, passes trough atomic positions in both layers 
(figure [T]), atoms in the 4 corners of the supercell are 
directly superimposed (figure [TJ. This corresponds to 
the so-called AA stacking. For (n, to = n + 1) cells, zone 
with AB stacking -graphite like stacking- are located at 
1 /3 and 2 /3 of the long diagonal (figure [T]) . 

Quite generally the structures that we study are fully 
characterized by the two index (n, m) of the rotation and 
also by the value of the translation vector between the 
two layers. However for sufficiently small values of the 
rotation angle 9 the Moire pattern depends essentially on 
the rotation angle and the translation of one layer only 
results in a rigid shift of the overall Moire structure. Also 
in the small angle limit, all dimensions of the Moire pat- 
tern scale like 1 /9 and in particular the size of locally AA 



or AB stacked regions and the distance between theses 
zones are proportional to 1 /9. 

Let us emphasize that the limit of small angle is dis- 
continuous: when the rotation angle decreases the dimen- 
sions of the Moire structure scale like 1/9 and therefore 
the structure does not tend to that for exactly zero an- 
gle. This singular, discontinuous, geometric evolution is 
associated to singular electronic properties as we show in 
this work. 

Ab-initio and tight-binding methods for electronic 
structure - Our approach combines ab initio calculations 
that can be used for unit cells containing up to 500-600 
atoms and tight-binding (TB) calculations that can be 
used for unit cells containing up to 15000 atoms or even 
more. Ab initio calculations are performed with the code 
VASP [23] and the generalized gradient approximation 
[25] . The C ultra soft pseudopotential [25] has been ex- 
tensively tested previously [TTJ, [2T] . The plane wave basis 
cutoff is equal to 211 eV. The empty space width is equal 
to 21 A, the interlayer distance is fixed to its experimen- 
tal value (~ 3.35 A) and no atom are allowed to relax for 
direct comparison with TB calculations. 

In the tight-binding scheme only p z orbitals are taken 
into account since we are interested in electronic states 
close to the Fermi level. Since the planes are rotated, 
neighbors are not on top of each other (as it is the case 
in the Bernal AB stacking). Interlayer interactions are 
then not restricted to ppir terms but some ppa terms have 
also to be introduced. For both terms, a parameter and 
a characteristic length have been fitted on AA, AB, (1,3) 
and (1,4) cells [27] to reproduce the ab initio dispersion 
curves according to: 

= -7 exp fq n \1 - ^j^j (3) 

Vp P<7 = 7i ex P f °* ( 1 - — ) J with — = — (4) 

a is the nearest neighbor distance within a layer a — 
1.42 A and ai is the interlayer distance a\ — 3.35 A. With 
this dependence, we take first neighbors interaction in a 
plane equal to 70 = 2.7 eV and second neighbors interac- 
tion in a plane equal to 0.1 x 70 [2S], which fixes value of 
the ratio q„/a. For undoped bilayers, all p z orbitals have 
the same on-site energy in both planes. We have cho- 
sen this on-site energy such that the energy of the Dirac 
point equals to for the results detailed here. 

Localization by Moire patterns - We have performed 
band structure calculations on a large number of struc- 
tures. For structures that could be studied by both ab 
initio and TB methods the agreement is always excellent. 
For example figure [2] presents the band structure calcu- 
lated with VASP and with the TB method for a small 
angle rotation: (n,m) = (6,7), 9 = 5.08° and N = 508. 
The two methods agree very well and show a decrease 
of the velocity by about 15 percent in that case. We 
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FIG. 2: (color on line) Band energy dispersion E(k) in 
commensurate bilayer cell (n, m) = (6, 7) for a rotation of 
= 5.08°. 
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FIG. 3: (color on line) Velocity ratio Vu /V morio for a com- 
mensurate (n, m) bilayer cell versus rotation angle 9: Circle 
VASP, cross TB calculations. Line is the model of Lopez dos 
Santos et al. [18]: V hi /V mon o = 1 - 9[i/(V mono K sin(f?/2))] 2 , 
with t = 0.11 eV and V mono K = 2j nV3 = 9.8 eV. 



have checked that the velocity remains isotropic in the 
linear region up to a few percents. For angles in the 
range 15° — 45° and for energies, of the order of a few 
meV (not shown here), the dispersion is not always lin- 
ear and a very small gap can even exist. Yet this small 
energy range is not experimentally relevant in most cir- 
cumstances. We note also that for angles within a few 
degrees of 0° or 60° the bands become flat and the energy 
range in which the bands are linear decreases. 

As stated above the geometry of the Moire depends 
essentially on the angle 9, in the limit of small 9. This 
reflects in the band dispersion which depends essentially 
on the rotation angle 9 between the two layers for all 
cases studied here and in particular in the small angle 
limit. We performed a systematic study of the renormal- 
ization of the velocity close to the Dirac point, compared 
to its value in a monolayer graphene, for rotation angles 
9 varying between 0° and 60° (figure |3| . 

The renormalization of the velocity varies symmetri- 
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FIG. 4: (color on line) Average participation ratio p(E) at 
energy E for several (n, m) bilayers. 



cally around 9 = 30°. Indeed, the two limit cases 9 = 0° 
-AA stacking- and 9 = 60° -AB stacking- are different, 
but Moire patterns when 9 — > 0° and when 9 — > 60° are 
similar because a simple translation by a vector trans- 
forms an A A zone to an AB zone. 

Focusing on angles smaller than 30° we define three 
regimes as a function of the rotation angle 9. For large 
9 (15° < 9 < 30°) the Fermi velocity is very close to 
that of graphene (see figure [3]) . For intermediate values 
of 9 (3° < 9 < 15°) the perturbative theory of Lopez dos 
Santos et al. |18j predicts correctly the velocity renormal- 
ization. But for the small rotation angles (9 < 3°) a new 
regime occurs where the velocity tends to zero. This re- 
markable localization regime cannot be described by the 
perturbative theory of Lopez dos Santos et al. |18j . 

In order to analyze this localization phenomenon we 
computed the participation ratio of each eigenstate 
defined by p(ip) = (N YV |(i|?/')| 4 ) _1 . \i) is the p z or- 
bital on atom i and TV is the number of atoms in a unit 
cell. For completely delocalized eigenstate p is equal to 
1 like in graphene. On the other hand, state localized 
on 1 atom have a small p value: p = l/N. The average 
participation ratio p at energy E is presented on figure 
H For intermediate values of 9 -bilayers (6, 7) in figure 
[4]- the participation ratio of states with energy close to 
is similar to that of other states of the bilayer. For very 
small 9 -bilayer (25,26) and (30,31) in figure [I]- states 
with energy around are strongly localized (the smallest 
p values down to 0.14). An analysis of spatial repartition 
of eigenstates, shows that theses states are localized on 
the AA zones of the Moire. This is illustrated on fig- 
ure [5J) showing 80% of the weight of an eigenstate with 
energy very close to in the (30, 31) bilayer cell. This 
new localization mechanism by the potential of the Moire 
leads to a peak in the local density of states of atoms in 
the AA zone (figure [5J)) whereas for Moire due to larger 
angle no strong localization occurs (figure [5^) . It would 
be interesting to check whether this localized peak can 
be observed in STS experiments. 
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FIG. 5: (color on line) TB Local density of states on atom 
located at the center of a AA type zone in (a) (6, 7) and (b) 
(30, 31) bilayers. Insert: Repartition of a electronic state at K 
with energy 0: Black small dots are the position of all atoms, 
red dots are atoms on which 80% of the electronic state is 
localized. 



Finally we have also studied the case of asymmet- 
ric bilayers with on-site energies that are different for 
atoms of the two planes. Such an asymmetry can occur 
in graphene multilayers due to the effect of doping. In 
the case of Bernal AB stacking an asymmetric bilaycr 
presents a gap but our TB calculations show that it is 
not the case for asymmetric rotated layers. In agreement 
with Lopez dos Santos et al. |18j we find that a difference 
between the on-site energies of the two layers results in a 
relative shift of the Dirac cones. The linear dispersion is 
conserved close to each Dirac point and no gap appears. 

Conclusions - In summary, by combining ab initio and 
tight- binding calculations, we have demonstrated that 
the dispersion relation stays linear in a rotated bilayer 
close to the Dirac point. Yet the velocity of electrons is 
rcnormalized. This renormalization depends essentially 
on a single parameter which is the rotation angle 9. For 
large 9, (15° < 9 < 30°), the Fermi velocity is very 
close to that of graphene. For intermediate values of 
9 (3° < 6 < 15°) the perturbative theory of Lopez dos 
Santos et al. |18j predicts correctly the velocity renor- 



malization. But for the small rotation angles (9 < 3°), 
a new regime occurs where the velocity tends to zero. 
This remarkable localization regime cannot be described 
by the perturbative theory of [TB]. In this regime elec- 
trons with energy close to Dirac point are localized in AA 
stacked regions. We believe that this localization regime 
should be observable since angles as small as a fraction 
of a degree occur in some rotated multilayers. 
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